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We explore the phase diagram of two-component bosons with Feshbach resonant pairing interac- 
tions in an optical lattice. It has been shown in previous work to exhibit a rich variety of phases and 
phase transitions, including a paradigmatic Ising quantum phase transition within the second Mott 
lobe. We discuss the evolution of the phase diagram with system parameters and relate this to the 
predictions of Landau theory. We extend our exact diagonalization studies of the one-dimensional 
bosonic Hamiltonian and confirm additional Ising critical exponents for the longitudinal and trans- 
verse magnetic susceptibilities within the second Mott lobe. The numerical results for the ground 
state energy and transverse magnetization are in good agreement with exact solutions of the Ising 
model in the thermodynamic limit. We also provide details of the low-energy spectrum, as well as 
density fluctuations and superfluid fractions in the grand canonical ensemble. 
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I. INTRODUCTION 

In the last few years there has been considerable ex- 
perimental and theoretical interest in studying Fesh- 
bach resonances between different atomic species and 
isotopes. This activity encompasses pairing interactions 
and molecule formation in a wide variety of Fermi-Fermi, 
Bose-Fermi, and Bose-Bose mixtures. Such systems pro- 
vide many possibilities for novel phases and phenomena, 
ranging from dipolar condensates [H Q to highly con- 
trollable chemical reactions Q . Recent examples include 
heteronuclear resonances in ^^Rb-^^Rb [j-H, ■*^K-®''Rb 
39K-87Rb [il, and 87Rb~"3Cs [ij Bose mixtures 
over a range of experimental parameters. 

Motivated by these developments, we recently investi- 
gated the phase diagram of two-component bosons pair- 
ing in an optical lattice [l^. Amongst our findings, we 
identified a paradigmatic Ising quantum phase transition 
occurring within the second Mott lobe. The principal aim 
of this manuscript is to provide a more detailed overview 
of this heteronuclear system, and to extend the scope of 
physical observables presented in Ref. We expand 

our previous exact diagonalization results in several di- 
rections, and confirm the additional Ising critical expo- 
nents, a = 0, 7 7/4 and 5 = 15. We also provide re- 
sults for the low-energy spectrum, as well as density fluc- 
tuations and superfluid fractions in the grand-canonical 
ensemble. We supplement this with a discussion of the 
evolution of the phase diagram with system parameters, 
and its direct connection to Landau theory. 

The layout of this manuscript is as follows. In Sec. [Tl] 
we describe the heteronuclear model with Feshbach in- 
teractions. In Sees, mil and IIVI wc discuss the mean field 
phase diagram and the Landau theory description. In 
Sec.|V]we focus on the Mott states and present a deriva- 
tion of the effective quantum Ising model. We confirm 
these findings in Sec. IVII bv exact diagonalization of the 
one-dimensional bosonic Hamiltonian. We conclude in 
Sec. lVIII We incorporate the principal results of Ref. [l5| . 



II. THE MODEL 

We consider two-component bosons with a "spin" in- 
dex ^, t, which may be different hyperfine states, isotopes 
or species. These components may form molecules, to, as 
described by the Hamiltonian 

iJ = ^ eaUia - X! X! (^Q^a + h.' 
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Here, Bose annihilation operators, where i labels 

the lattice sites, = a\^aia, and a =i,\, m; Eq, are on- 
site potentials, ta are hopping parameters, (ij) denotes 
summation over nearest neighbor bonds, and Uaa' are 
interactions. Molecule formation is described by the s- 
wave interspecies Feshbach resonance term 



h.c). 



(2) 



Similar problems have been studied in the continuum 
limit with s-wave [l^ and p-wave resonances • Closely 
related homonuclear systems have also been considered in 
the continuum 18 22j and on the lattice [23l426| . Normal 
ordering implies : niaUia '. = nia{nia — 1) for like species, 
and : riiaTiia' ■ = Uiariia' for distinct species. For sim- 
plicity we consider hardcore atoms and molecules and 
set J7mt = Umi = U and U^i = V. We begin work in the 
grand canonical ensemble with = H — ii^^Nt^ — ^dNd, 



where iVr 



-2nj„i) is the total atom num- 



ber, including a factor of two for molecules, and A^d = 
^^(n,;^ — n.4) is the up-down population imbalance. The 
chemical potentials may be absorbed into the coefficients, 
e-|- = ei--MT-MD, ej, = ei-^T + MD, and e„, = e„-2/iT- 
For the numerical investigation of the Ising transition we 
shall subsequently switch to the canonical ensemble with 
a fixed total density, = + + 2nm = 2. 
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FIG. 1. Zero hopping phase diagram showing the minimum 
energy eigenstates in tlie basis \ni,nf;nm). The Feshbacli 
coupling mixes |1, 1; 0) and |0, 0; 1) to yield |±) with energies 
E±^im- h/2 ± VW2F+52, where h = - €^ - e-f - V . 
The other states have E{ni, 

Unm{n-f + n_i). The total density pT = ni + n-f + 2nm is pinned 
to integer values and increases with fiT- Increasing (decreas- 
ing) /iD favors up (down) atoms. The topology changes with 
the system parameters and depends on the signs of the energy 
gaps, A±, defined in the text. We set ej, = e-f = g = 1, [/ = 
and (a) e„ = 2, F = 0, (b) = 3.5, V ^ 1, (c) e„ = 3, 
V = 1.2, (d) em = 2.5 V = 1. The panels shown in Fig. [2] 
correspond to vertical slices through diagram (a). 



III. PHASE DIAGRAM 

As discussed in Ref. [l^, in elucidating the zero tem- 
perature phase diagram the zero hopping limit provides 
a useful anchor point. In particular, the topology of the 
zero hopping phase diagram depends on the Hamiltonian 
parameters via the energy difference between proximate 
phases; see Fig.[TJ In this respect it is convenient to define 
the energy gap, A_ = E{0, 0; 0) — E- , between the vac- 
uum state and the second Mott lobe, evaluated at the in- 
tersection point E{1,0; 0) = E{0, 1; 0) = £^(0, 0; 0), where 
we label the energies E(n^, n-^; n„i) in the occupation ba- 
sis. Similarly, we define A+ = E^ — i?(l, 1; 1), between 
the second and upper Mott lobes, where i?(l,0;l) = 
E{0, 1; !) = £;_. This yields A± = y^g^ + {h/2)^ - V± 
h/2, where ft, = e,„ — ej^ — e-^ — y. One obtains the struc- 
ture shown in Fig. [ija) for A± > 0, (b) for A_|_ > and 
A_ < 0, (c) for A± < 0, and (d) for A+ < and A_ > 0. 

For simplicity, we begin with the parameters used in 
Fig. [Ija), where U ^ V = 0. Since we include Feshbach 
resonant interactions, g, this limit captures many of the 
principal features and phases of the interacting problem, 
including the presence of Mott states. We will incor- 
porate the effects of finite U and V in the subsequent 
discussion. The mean field phase diagram is obtained by 
minimizing the effective Hamiltonian 
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FIG. 2. Evolution of the mean field phase diagram with /iD, 
corresponding to vertical scans through Fig. [TJa). We set 
€^ = ei = \, em = 2, g = I, U = V = = ti = t, 

and tm = t/2. We indicate the one-component up, down 
and molecular superfluids, by SF-f, SFj, and SFm, while MI 
denotes a Mott insulator. Phase 3SF has all three components 
superfluid. We denote first order (continuous) transitions by 
double (single) lines. Junctions between phases are indicated 
by a dot, and the termination of first order lines by a cross. 



where is the single site zero hopping contribution to 
([T|), z is the coordination, and (t>a = {a-ia)', see Fig. [51 
The phase diagram is symmetric under /ix — > 2 — due 
to invariance of the Hamiltonian ([1]) under particle-hole 
and spin flip operations, aa ^ al^, iit ^ + U — ht, 



a-|-, when ~ t-f and h ~ 0; this extends the 



"top-bottom" symmetry in Fig. [TJa). Likewise, the sys- 
tem is invariant under /id ^1' ~ ~ Md, a-f 
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FIG. 3. Mean field phase diagram in the presence of inter- 
actions, U and V . We set e-|- = = 1, £,„ = 2, g = 1, 
M — ^ir~ ^5 = and /id = 0.45 as in Fig. [2jc), and con- 
sider the efi^ects of U and V in turn, (a) [/ = 1, = 0, show- 
ing the increase in extent of the |— ) Mott lobe, (b) [/ = 0, 

V — 1, showing the appearance of asymmetry, (c) U = 1, 

V — 1 showing the combined effect of both interactions. The 
characteristic phases and phase transitions mirror those seen 
in the U = V = limit shown in Fig. (2] 



and ti t^, which extends the "left-right" symmetry in 
Fig. [TJ The phase diagram has a rich structure and ex- 
hibits distinct Mott insulators, single component atomic 
and molecular condensates, and a phase with all three 
species supcrfluid. However, two-component supcrfluids 
are absent due to the structure of the Feshbach term 
The main panel shown in Fig.[21[c) also displays 
an intricate network of continuous quantum phase transi- 
tions, quantum critical points and first order phase tran- 
sitions, where the latter are inferred by discontinuities in 
the order parameters and derivatives of the ground state 
energy. In particular, the first order segments shroud the 
second Mott lobe and overextcnd beyond the junctions 



of the proximate supcrfluid phases. Similar features also 
emerge in the two-component Bose-Hubbard model in 
the absence of Feshbach interactions [l^. As may be 
seen by tracking the evolution with /.iD, these first order 
segments emerge from an underlying tetracritical point, 
as shown in panels (a) and (b). In a similar way, the 
tetracritical points shown in Fig. [H^c) may bifurcate into 
segments connected by first order transitions as shown in 
panels (d) and (e). As we shall discuss in Sec. |lVl the 
presence and transmutation of these elementary critical 
points and first order segments may be seen from a re- 
duced two-component Landau theory. Before embarking 
on this discussion, let us note that these principal fea- 
tures also emerge for non- vanishing U and V , as shown 
in Fig. [31 Although finite interactions induce quantitative 
distortions of the phase diagram depicted in Fig. [21 the 
characteristic phases and phase transitions arc nonethe- 
less present in the U = V = limit. 



IV. LANDAU THEORY 

To understand the structure of the mean field phase 
diagram we develop a Landau theory description. The 
eigenenergies of the effective single site Hamiltonian ([3]) 
may be calculated by treating the off-diagonal hopping 
contribution as a perturbation, 



V 



E 



(4) 



where a =4,, t, ^n. In principle, this may be performed to 
arbitrary order using the general formalism in Ref. [28j . 
To illustrate the observed topology it is sufficient to ob- 
tain the appropriate eigcnencrgy to fourth order 



Er = E^ + {r\VS^V\r) 
+ {r\VSlVSlVSlV\r) 



{r\VSi.VSi.V\r) 
{r\VS^V\r){r\VS^V\r} 



(5) 



where r labels the zero hopping Mott state, and the op- 
erator = Y^n^r \n){n\l{E^r - is introduced [H]. 
The expansion ([5]) takes on a simplified form due to the 
off-diagonal nature of the perturbation. One obtains 
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+ \Y.>^o.p\<t>^?\<t>p?+o{<t>^), 



h<t>i + li-c.) 



(6) 



a/3 



where the explicit parameters, but not the overall struc- 
ture, depend on the unperturbed Mott state. The tran- 
sition to the one-component superfiuids is determined by 
the vanishing of the quadratic mass terms; expressions 
for these coefficients are given in App. [XI The function 
([6]) is minimized when aT:g(^4'*^4>^(j)^) = tt, and with- 
out loss of generality we may take the fields to be real. 
To see how a Landau theory of this type gives rise to 
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FIG. 4. Topologies of the reduced Landau theory ([8]) with 
a, 13 = 1,2 and a term {4>i + (fit)/^ added for stabiUty. We set 
All = A22 = 1 and consider the evolution with A12. We 
denote first order (continuous) transitions by double (sin- 
gle) lines. Junctions between phases are indicated by a 
dot, and the termination of first order lines by a cross. (I) 
A fa < A1 1A22: tetracritical point for A12 = 0.8. (II) A12 > 
\/AiiA22: first order transition between one-component con- 
densates for A12 = 1.5. (Ill) A12 < — \/AiiA22: first or- 
der transition with condensation of both order parameters 
for A12 = -1.5. 



the observed evolution of the tetracritical points we con- 
sider the situation where one of the fields is massive with 
m » 0. Replacing this field by its saddle point solution 
one obtains the energy in terms of the two remaining 
fields. As a concrete example let us consider the vicinity 
of the lower tetracritical point shown in Fig.[2l^c), where 



> 0. 
given by 



The value of ipi which minimizes the energy is 



dE 



(7) 



Taking (f)^ = ~l4'm<t>'[/'nii + 0(0"*) we obtain a reduced 
two-component Landau theory 



E = E. 



0{<P% (8) 



Q/3 



where a, /3 =t, rn, and the interspecies density-density in- 
teraction has been renormalized to K^m = '^f™ ^ 7^/™4. 
whilst Aqq. = Aqq remains unchanged. The behavior of 
this Landau theory is governed by the sign and magni- 



tude of A 



There are three distinct cases to con 
A, 



> A^„j this 

vara 



sider as illustrated in Fig. |4l (I) if K^^i^mm 
describes a tetracritical point; (II) if A-j-^ > 
then there is a first order transition between the two sin- 
gle component superfiuids; (III) if K^^a < — \J A^^^mm 
then the system is unstable at fourth order indicating 
the presence of a first order transition to condensation 
of both order parameters. For the parameters used in 
Fig. [2] we find explicit examples of types (I) and (III); 
see Figs. [5] and [6l At the lower tetracritical point in 
Fig. Wis), where = m™ = 0, the Landau coefficients 
are given by A-t-|- = 4e-|-, 
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FIG. 5. Magnified portion of the lower left region of Fig. [2] 
showing the underlying tetracritical points as /id is varied. 
The lower and middle tetracritical points bifurcate at fiu ~ 
0.203 and /id ~ 0.445 respectively. The former corresponds 
to the reduced Landau theory criterion A^^ — — ^ A-|-|-Am»n 
obtained by saddle point elimination of in the perturbation 
expansion around the vacuum state, |0,0;0). The latter cor- 
responds to A^m ~ — i/A^^Amm after eliminating in the 
perturbation expansion around the second Mott lobe, |— ). 
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FIG. 6. Magnified portion of the second lobe tip in Fig. [2] 
showing the retreat of the first order transitions and the emer- 
gence of a tetracritical point at /id ~ 0.715. This is given by 
the reduced Landau theory criterion, A^^ = — ^ A^^ A_l 



ob- 



tained after saddle point elimination of 
theory around the second Mott lobe, |— ). 



in perturbation 



(9) 
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and 



A-|-„i — — 



+ - 



(10) 

where = (e^ + et)/^ + — \/'ig^ + J^^)/3 for 2tm = 
'•i ~ ~ ^' = h + {l^ — e^) /4 and £ = + e,„. Setting 



A- 



as appropriate for Fig. [2J the bifurcation condition 

0.2035 + (e^ 



mm 



yields 



^i)/2. 



This is in quantitative agreement with the bifurcation 
point shown in Figs. [2] and EJ where 5 = 1. In a similar 
fashion, for the parameters used in Fig. [2l the |— ) lobe 
tetracritical points bifurcate at /^d ~ 0.445 and /iD = 
(4/v^- 1)1/2 ^ 0.715; see Figs. [5] and [H 



V. MAGNETIC DESCRIPTION 

Having discussed the phase diagram we turn our atten- 
tion to the Mott state with total density = 2. This 
reveals Ising transitions in both the heteronuclear and 
homonuclear lattice problems (l5| . To explore this sec- 
ond Mott lobe, where the Feshbach term is operative, 
we adopt a magnetic description. With a pair of atoms 
or a molecule at each site, we introduce effective spins 
|4J.) = |1,1;0), Itr) = |0,0;1); see Fig. [H The operators 



S+ 



1 



(11) 

or equally 5^ = [n,n — [n^ + 72-|-)/2]/2, form a representa- 
tion of su(2) on this reduced Hilbert space. Although the 
spin representation (jlip contains three bosons it is cl osely 
related to the more familiar Schwinger construction [30[ . 
Deep within the Mott phase we perform a strong coupling 
t/U expansion [3l|: 

H = JY.sts^^+Y. (^^^ + r^f ) ' (12) 

where F = 2g, ft, = — ej, — e-f- — and the exchange 
interaction is given by 
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FIG. 7. Schematic representation of the second Mott lobe 
with pT = 2 showing the presence of a pair of atoms |L or a 
molecule 1^ at each site. XY exchange involves interchanging 
a molecule and a pair of atoms, and is suppressed relative to 
the Ising interaction for small hoppings. 



the form of a quantum Ising model in a longitudinal and 
transverse field. The longitudinal field, h, refiects the en- 
ergetic asymmetry between a molecule Iff) , and a pair of 
atoms 14)- The transverse field, F = 2g, encodes Fesh- 
bach conversion and induces quantum fluctuations in the 
ground state; see Fig. [T] In particular, in one dimension 
and with h = 0, the model ([T^ is exactly solvable by 
fermionization. It exhibits a quantum phase transition 
at F = J/2 from an ordered to disordered phase js^-fsj]- 
For h 0, the location of this transition is modified; see 
for example Ref. [ssj . 



VI. NUMERICAL SIMULATIONS 

The model (|12[) plays an important role in quantum 
magnetism and quantum phase transitions [sgI. |37{. To 
verify this realization in our bosonic model ((!]), we per- 
form exact diagonalization on the ID quantum system ([T]) 
with periodic boundary conditions, at zero temperature. 
The large Hilbert space (x 2"^^ of the multicomponent 
system under consideration restricts our simulations to 
N < 8 sites. In Sec. IVI Al we present results for the 
overall phase diagram in the grand canonical ensemble 
(allowing all possible states) for N = 6 sites, before mov- 
ing on to a detailed canonical ensemble study (allowing 
only states with fixed total density pT = 2) of the Ising 
quantum phase transition in Sec. IVI Bl Comparing di- 
rectly to results for the Ising model at a given system 
size provides a good understanding of finite-size effects 
in the bosonic problem. 
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t| + i| 



U-V 2U 



(13) 



Here we focus on the antiferromagnetic case with J > 
0. In writing (|12p we have omitted the constant, — 
h/2 — Jz/8 per site which is necessary for quantitative 
comparisons; see App. [BJ More generally, in the presence 
of asymmetry between U„i^ and U,ni 



^1 



7nX 



U^i Umf — Ufl Umi + U, 



mt 



(14) 

The structure of this exchange is readily seen from the 
energy cost for each individual hopping process in second 
order perturbation theory. The Hamiltonian takes 



A. Grand Canonical Phase Diagram 



In Fig. [5] we prese nt results for t he fluctuations in the 
total density, a = -y/ (n^) — (nT)^ where nx = J2a "-i" 
for an arbitrary site i, and the total superfluid fraction 



N 



92 



(15) 



Here we impose a phase twist & <^ n and calculate 
the change in the ground state energy i5gs(0) [H, [soj . 
These key observables show the onset of superfluidity, 
and the results are in good qualitative agreement with 
each other, and with the mean field phase diagram 
shown in Fig. [21a,). As found in the single component 
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FIG. 8. (Color online) (a) Density fluctuations and (b) super- 
fluid fraction for the ID bosonic model ([T]) obtained by exact 
diagonalization for N = 6 sites. We set ej, = e-|- = 1, = 2, 
U = 1, V = 0, t = = tf = 2t,n and /id = 0.45. The extent 
of the Mott lobes is in good agreement with the mean held 
predictions in Fig. EJa). The system sizes are insufflcient to 
resolve the quantum phase transitions between the distinct 
superfluids shown in Figs. [2] and [S] 



Bose-Hubbard model [40| , the Mott lobes develop sharp- 
ened tips due to enhanced fluctuations in one dimension 
[111, 11^. However here, the Mott lobe remains symmetric 
around the mid-point due to the symmetry of the Hamil- 
tonian ([Ij under particle-hole and spin flip operations, 

ac O- ajj, Cm + U — ht, aj. ^ oti ^'^^ '^^^ chosen 

parameters with = t-^ and h = 0. The system sizes 
accessible by exact diagonalization arc insufficient to re- 
solve the quantum phase transitions between the distinct 
superfluids shown in Figs. [2] and [31 However, as we will 
demonstrate in Sec. IVI Bl they provide compelling evi- 
dence for the magnetic structure of the Mott phases [l^ . 



B. Ising Transition in the Second Mott Lobe 

To employ greater system sizes we switch to the canon- 
ical ensemble with fixed density, px = 2. To study the 
Ising transition we work in the region of small hopping 
parameters, and begin with h = before exploring fi- 
nite fields. Due to the absence of spontaneous sym- 
metry breaking in finite systems, the staggered mag- 
netization vanishes in the absence of an applied stag- 
gered field. In view of this we focused our previous nu- 
merical investigation [l^ on the pseudo staggered mag- 
netization, m = {\T,ii-'^yS-\)/N m, where S'f = 
[nim — {ni^ + nii)/2\ /2 and additional modulus signs are 
incorporated. As shown in Fig. IHl^a), this quantity is 
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FIG. 9. (a) Pseudo staggered magnetization, m, versus F = 
2g, obtained from exact diagonalization of the ID bosonic 
model lU for different system sizes. A'', at density pT = 2. We 
take e_i = — 1, em = 2 (corresponding to h = 0), U = 1, 
V = and set t = ti = tf = 2tm = 0.01 corresponding to 
J = nt^ /A. (b) Finite size rescaling of mN^ versus T = 2g 
showing the presence of the Ising quantum phase transition at 
r = Fc « 2.08 X lO""* and the critical exponent /3 = 1/8 of the 
2D classical Ising model. The critical coupling Fc ~ 0.49J is 
close to the small hopping Ising result, Fc = J/2, (c) Scaling 
collapse as a function of (F — Fc)A^'''^ corresponding Xo v = 1. 



rendered finite. Adopting the finite-size scaling form, 
m = N-P/^rh [(r - Tc)N^/''] [H, we plot mN^/^ ver- 
sus g for different system sizes, A^, in Fig. [Hlb). The 
curves cross close to the critical coupling, Fc = J/2, of 
the purely transverse field Ising model. Moreover, the 
scaling collapse shown in Fig. is consistent with the 
critical exponents /3 = 1/8 and = 1 for the 2D classical 
model. In spite of this evident success, it is clearly desir- 
able to examine this transition in direct physical observ- 
ables, and we turn our attention to this below. In par- 
ticular, we establish a direct connection with analytical 
results and confirm additional Ising critical exponents. 

As shown in Fig. [TUl the correlation length exponent, 
V ^ 1, also follows directly from the gap data. Here we 
plot {A/J)N'' versus F = 2g, where A = E^^ - Eg^ is 
the energy gap between the first excited state and the 
ground state. As indicated in Fig. [TUTd). the values of 
the gap are also in very good quantitative agreement with 
finite-size simulations of the Ising model. 

In a similar fashion we obtain the exponent, 7 = 7/4, 
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FIG. 10. (a) Dimensionless energy gap A/ J versus T = 2g 
obtained by exact diagonalization of the ID bosonic Hamil- 
tonian for different system sizes, A^. We set ei — e-f — 1, 
= 2, U = 1, V = Q, t = ti ^ = 2tm ^ 0.01. (b) 
Scaled energy gap {A/J)N" versus T = 2g sfiowing tfie pres- 
ence of the Ising quantum phase transition at the numeri- 
cally extracted value, Fc ~ 2.08 x 10"*, and the correlation 
length exponent v = 1. (c) Scaling collapse as a function of 
(F — Fc)A'^^'^'' corresponding to v = I. (d) The numerically ex- 
tracted critical coupling, Fc « 2.08 x 10"'*, obtained from the 
crossing point in Fig. [TD](b) is consistent with the gap closing 
in the thermodynamic limit. A'' — >■ cxj. We set F equal to mul- 
tiples of J/2 (or the numerically extracted bosonic coupling, 
Fc) for the Ising (bosonic) model. 



from the magnetic susceptibility x — c^"^st/9/ist|/i_>o+ ' 
where mgt = {J2ii~^y^i) ^^'^ ^st is a staggered 
magnetic field applied to the bosonic Hamiltonian ([1}, 
AH = -hstE^i-'^yS!■, see Figs. Efa) and (b). A sig- 
nificant advantage of this somewhat more involved proce- 
dure is that the staggered magnetic field couples directly 
to the genuine order parameter, mst, without the need for 
modification. Likewise, we compute the transverse sus- 
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FIG. 11. (a) Staggered magnetic susceptibility obtained by 
exact diagonalization of the ID bosonic Hamiltonian IT} re- 
vealing the Ising quantum phase transition at Fc ~ 2.08 x 
10^^, and the critical exponent 7 = 7/4. We set e_i = €^ — 1, 
e,r, = 2, U ^ 1, V ^ 0, t ^ ti = = 2tm = 0.01. (b) Scaling 
collapse showing u — 1. (c) The transverse magnetic suscepti- 
bility, xx, plays an analogous role to the specific heat capacity 
of the 2D classical Ising model, (d) The logarithmic diver- 
gence at the critical point yields a = (e) log(mst) ver- 
sus log(/ist/J) to extract the exponent S = 15 for the bosonic 
model Q using m^t ~ h^^ for F = Fc ~ 2.08 x 10"-*. 
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J 

FIG. 12. Ising quantum phase transition in the bosonic 
model (HI showing the emergence of the critical line, Fc = J/2, 
in the regime of small hoppings. We set £4, = e-j- = 1, em = 2, 
U — 1, V — 0, t — ti — tf = 2tm and vary the hopping, t, 
corresponding to different values of J = 17t'^/4. 



ceptibility, x± = dm±/dT, where m± = J2i{^i)/^ is 
the transverse magnetization. Since the transverse field 
acts to disorder the system, this plays a similar role to 
the specific heat capacity of the 2D classical Ising model 
[43} . The dependence shown in Fig. HHd) is consistent 
with the Ising critical exponent a = 0. At criticality, the 

relation, mgt ^ht^ ^ ^^^'^ furnished with (5 = 15 as 
shown in Fig. [Tlje). 

All of the above diagnostics confirm the consistency 
of the measured Ising critical point [isj . and we may 
track this transition within the Mott lobe. The results 
in Fig. [12] show clear Ising behavior at small hoppings 
corresponding to Fc ~ J/2. With increasing hopping the 
boundary peels away from this linear slope and we find 
Fc < J/2. Additionally, enhanced finite-size effects result 
in slightly different estimates of Fc from to, A and x; see 
Fig. 1121 This will be explored in future density matrix 
rcnormalization group (DMRG) work p^ . 

As shown in Fig.[T3](a), the results obtained for = 8 
sites are in good agreement with exact results pertaining 
to the thermodynamic limit of the quantum Ising model 
for both the ground-state energy density [U 



/ (ifcV4r2 

47r J a 



J2 +4FJcosfc, 



43 



and the transverse magnetization [3. 

dk 2F + Jcosfc 



(5r> = - 



27r V4F2 + J2 +4FJcosfc 



(16) 



(17) 



In this comparison we define E' = Egg — (e,„ — h/2 — 
Jz/8)N for the bosonic system |T|) in order to take into 
account the constant offset in the mapping (|12p. The 
residual finite size effects may be used to find the central 
charge of the bosonic system. With periodic boundary 
conditions the ground state energy depends on the sys- 
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FIG. 13. Comparison of (a) the shifted ground state energy 
Egg = Egs — {Sm — h/2 — Jz/8)N, and (Inset) transverse 
magnetization m± = Xli('S'f)/N of the bosonic model ([T]) 
and the exact thermodynamic results (|16p and (|17p for the 
Ising model. We set = = 1, em ~ 2, U — 1, V = 0, 
t = ti = — 2tm = 0.01. (b) Finite size corrections for iJgg 
of the bosonic model {TJ, with F = Fc ~ 2.08 x 10"''. We use 
the finite size scaling result, E'^^/JN = e'^ / J — ■kc{v / J) / QN"^ . 
The intercept at -0.3164 is in good agreement with the exact 
Ising result (|16|l in the thermodynamic limit, e'^/J — — l/vr « 
-0.3183 when F = J/2. The slope, -7rc/12, yields the cen- 
tral charge, c « 0.52, where we use the characteristic veloc- 
ity V — J/2, (c) Comparison of the low-energy spectra ob- 
tained by exact diagonalization for N = S sites in the bosonic 
model ((!]) with the same parameters as used in panel (a) with 
T = 2g = 2J, and the Ising model gl]) with J = 17f^/4, 
ft = 0, and F = 2J. 



thermodynamic limit, v is the effective velocity in the lin- 
earized low energy dispersion, and c is the central charge 
[46l-l4i|. For the ID Ising model the dispersion re- 
lation may be obtained by fermionization |32l434| . This 
yields the characteristic velocity, v ~ dek/dklf^^^ = J/2. 
The recovery of the Ising model central charge, c = 1/2 
[HjH^I, is shown in Fig. [T51' bV Although the majority of 
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FIG. 14. (a) Evolution of the Ising quantum phase transi- 
tion with the longitudinal field, h, in the ID bosonic model 
([T]). Here = 2 + /i and all other parameters as in Fig. [9] 
(b) Quadratic correction to the ground state energy for the 
bosonic Hamiltonian ([l} due to a longitudinal magnetic field, 
h, in the vicinity of the h = Ising critical point. We set 
e:^^e^ = l,e^ = 2,U^l,V = 0,t = tl = ^ = 2tm= 0.01 
and r = Fc ~ 2.08 x 10"". The results are close to the theo- 



retical value, SEi 



(2) 



-0.0706h^N, at F = Fc [SE 



our discussion has been on the ground state properties, 
and the first excited state, the agreement between the 
Ising model description and the bosonic Hamihonian also 
extends to the low-energy spectrum shown in Fig. II3f c). 

Having provided a thorough description of the Ising 
quantum phase transition in the absence of a magnetic 
field, let us now consider its effect. In the presence of a 
finite longitudinal field, h, the location of the Ising criti- 
cal point is modified as shown in Fig. fW a). The results 
obtained from the pseudo staggered magnetization, mag- 
netic susceptibility, and gap data are in good agreement 
with exact diagonalization and DMRG results [1^ for 
the Ising model. Finite-size effects increase with increas- 
ing magnetic field, and for h/J > 0.6 and = 4,6,8, 
clean scaling is observed for m, but not for A and x- 
In addition, departures from the h = critical point 
yield the expected quadratic dependence of the ground 
state energy on the longitudinal field, h. The relation 
SE^V = -OmOGh^N [3l,|5l| is recovered in Fig.Mi>)- 
These results show the persistence of an Ising quantum 
phase transition in the bosonic model without the need 
to fine tune to zero magnetic field. 



VII. CONCLUSIONS 

In this work we have provided a detailed study of 
bosonic heteronuclear mixtures with Feshbach resonant 
pairing in optical lattices. The model displays a rich 
phase diagram with an intricate network of quantum 
critical points and phase transitions, and we have an- 
chored this behavior to the predictions of a reduced two- 
component Landau theory. We have substantiated our 
previous findings of an Ising quantum phase transition 
occurring within the second Mott lobe by signifi- 
cantly extending the range of physical observables. In 
particular, we have confirmed the additional Ising crit- 
ical exponents, a = 0, 7 ~ 7/4, and i5 = 15, for the 
one-dimensional bosonic system. There are many direc- 
tions for further investigation including the superfluid 
properties, and studies away from the hardcore limit. 
Cold atoms in optical lattices may represent ideal sys- 
tems in which to explore magnetization distributions [52 1 
and quantum quenches [53|.[m| in quantum Ising models. 
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Appendix A: Landau Coefflcients 

To calculate the phase boundaries delimiting the single 
component superfluids from the Mott lobes it is sufficient 
to examine where the relevant mass terms change sign. 
In particular, for the second Mott lobe, |— ), we find 



i5_ =i?^i + i^m,02 



(Al) 



where 



h(h-e -U) 



X [(X + Uy - (I - 




x{x-P) + t^{eu + a) 



Here we denote x = Vg^^+JhJW. a^h/2+ {U + V)/2 
and (i^{V - U)/2. The coefficients for |0, 0; 0) read 



mt^(0,0;0) = ztt4 ( 1 - ^ 



m™(0, 0;0) = ztmil 



(A4) 
(A5) 
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whilst those pertaining to |1, 0; 0) read 

1 _ ^^'i 



to™(1,Q;0) = zt„ 
m;(l,0;0) = zi; 

mt(l,0;0) = zi-f I 1 



u 



h + e^ + V 



x'~ii + ~et + vr 



(A6) 
(A7) 

(A8) 



The remaining coefficients may be obtained by inter- 
changing particles and holes and ups and downs as appro- 
priate. For example, the coefficients TOq,(0, 1;0) are ob- 
tained by interchanging io-f in equations (|A6p — (|A8p . 
Likewise, the coefficients mQ,(l, 1; 1) are obtained by the 
particle-hole transformation, e-^i — >■ — (e^j. + U + V) and 
e„i -(e,„ + 2U), on equations (|A4|) and (|A5|) . 



The only process contributing at second order is a parti- 
cle hopping to a neighboring site and back: 



am (B6) 

Noting that 

aL«j™^o = |l,l;l).|0,0;0),(^j4.|7'o, (B7) 

a\^a,^Va = |0,l;l),|l,0;0),(%t.|^o, (B8) 

aW^o = |l,0;l).|0,l;0>,(4,^.|Po, (B9) 

one obtains 



Appendix B: Derivation of the Ising Hamiltonian 

In this Appendix we derive the effective spin model de- 
scribing the second Mott lobe of the heteronuclear Hamil- 
tonian. Employing the spin operators as defined in equa- 
tion pip and noting the form of the zero hopping eigen- 
states we may re-write the bosonic Hamiltonian ([T|) as 



H = NEoVo + H" + H', 
where N is the number of lattice sites. 



and 



E 



EUaa' 
9 



(Bl) 
(B2) 



H'= '^xk.X^d- E Ual^ja+h.c). (B3) 



J,fT — ± 



Here, Eq = ern-h/2, and T^o = 11, \h){h\ + \1l^){1\i\ is the 
projection operator onto the quasi-degenerate subspace 
of a molecule or two atoms on each site. The effective 
Hamiltonian within this subspace is [28| 

Hs = NEoVo+VoH'Va+VoH' S'^ H'Vo+0{H'^), (B4) 

where S'^ = J2m^t,iy \m){m\/{NEo - Since only 

the hopping contribution to H' connects this subspace to 
the other states one immediately finds 



i,cr—:L 



(B5) 



VoH'S' H'Vo = 



2U V-U 



(BIO) 



where Um-t ~ Umi ~ U and U^^ ~ V. Using the identity 



IM.>(Md = (l/2 + 5|)(l/2-5f) 



(Bll) 



one obtains 



VoH'S'' H'Vo = ^VoJj2 (j - ^f^l) ^0, (B12) 



where J is given by equation ([T^. We conclude that 



H^Vn 



(B13) 



Vo, 



as stated in equation p2|) . In the zero hopping limit 
the eigenstates of Hamiltonians (HI) and (|B13p coincide, 
including the mixing induced by the Feshbach coupling. 
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